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Abstract 

Two-field slow-roll inflation is the most conservative modification of a single-field model. The 
main motivations to study it are its entropic mode and non-Gaussianity. Several years ago, for 
a two-field model with additive separable potentials, Vernizzi and Wands invented an analytic 
method to estimate its non-Gaussianities. Later on, Choi et al. applied this method to the model 
with multiplicative separable potentials. In this note, we design a larger class of models whose non- 
Gaussianity can be estimated by the same method. Under some simplistic assumptions, roughly 
these models are unlikely able to generate a large non-Gaussianity. We look over some specific 
models of this class by scanning the full parameter space, but still no large non-Gaussianity appears 
in the slow-roll region. These models and scanning techniques would be useful for future model 
hunt if observational evidence shows up for two-field inflation. 
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I. INTRODUCTION 



Cosmic inflation [1] is a great idea to solve some cosmological problems and to predict the 
fine fluctuations of cosmic microwave background (CMB). Hitherto the surviving and most 
economical model of inflation involves a single scalar field slowly rolling down its effective 
potential 0, Isj , with a canonical kinetic term and minimally coupled to the Einstein gravity. 
We will call it the simplest single-field inflation, although there is still freedom to design its 
exact potential. The single-field inflation passed the latest observational test [4J successfully, 
even with the simplest quadratic potential. 

Nevertheless there are perpetual attempts to modify the simplest single-field inflation. 
Some of them are motivated by incorporating inflation model into certain theoretical frame- 
works, such as the standard model of particle physics [!, @] or string theory Q. Some others 
put their stake on signatures that cannot appear in the simplest single-field model, such as 
a large deviation from the Gaussian distribution in the CMB temperature fluctuations. 1 

Among these modifications, the two-field slow-roll inflation is the most conservative one, 
at least in my personal point of view. It introduces another scalar field rather than a non- 
conventional Lagrangian such as non-canonical kinetic terms or modifications of gravity. 
It also retains the slow-roll condition, which makes the model simple and consistent with 
the observed CMB power spectrum. If both conventional Lagrangian and non- conventional 
Lagrangian are adaptable to the observational data, then the model with conventional La- 
grangian would be more acceptable, unless there are better and solid theoretical motivations 
for non-conventional Lagrangian. 

On the observational side, two new features arise in two-field model. First, the model is 
able to leave a residual entropic perturbation between the fluctuations of dark matter and 



CM B |21| . |22j. Second, in a simple model with quadratic potential, numerical computations 
23, 24 found that the non-Gaussianity can be temporarily large at the turn of inflation 
trajectory in field space. Longer-lived large non-Gaussianities were discovered recently by 
[25H27I] in many other two-field models. 2 



Compared with the simplest one-field inflation, the field space becomes two-dimensional 
in a two-field model. When the inflation trajectory is curved in field space, the entropic 
perturbation will be coupled to the adiabatic perturbation. So there are more uncertainties 
in calculation of cosmological observables, such as power spectra of CMB and their indices. 
It would be more complicated to honestly compute the bispectra and non-linear parameters, 
which reflect the non-Gaussianity of the primordial fluctuations. 



Fortunately, based on the extended 8N- formalism [33|, Vernizzi and Wands [24J invented 
an analytic method to estimate such non-Gaussianities. They demonstrated the power of 
this method in a two-field model with additive separable potentials. This method was later 



applied by Choi et al. |34| to a model with multiplicative separable potentials. 

Encouraged by the method of Vernizzi and Wands, we tried to improve it for the two-field 
slow-roll model with generic potentials but failed. Finally, we only designed a larger class 



1 As a partial list, see [8M20I] and references therein for various models and recent development along this 
direction. 



The readers may refer to [28] for a review on this topic, and to [29, 30 [ for pioneer works that computed 
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analytically the non-Gaussianity expected in multi-field inflation. By studying the loop corrections, 
obtained an observable level of non-Gaussianities, even when the two-field model is of the slow-roll variety 
with canonical kinetic terms and in the framework of Einstein gravity. 
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of models whose non-Gaussianity can be estimated by this method. It is a class of models 
whose potential take the form W(w) with w = U((p) + V(x) or w = U(ip)V(x)- Here W(w), 
U((f) and V(x) are arbitrary functions of the indicated variables as long as the slow-roll 
condition is satisfied. Scalar fields ip and x are inflatons. 

The outline of this paper is as follows. In our convention of notations, we will prepare 
some well-known but necessary knowledge in section [TT] concisely. In section IHH we will 
present the exact form of our models, whose non-linear parameters will be worked out in 
sections IIVI and [V] Some specific examples are investigated in section IVI1 We summarize 
the main results of this paper in the final section. 



This is a note concerning references [24], |34J . Some of our techniques stem from these 



references or slightly generalize theirs. Sometimes we employ the techniques with few expla- 
nation if the mathematical development is smooth. To better understand them, the readers 
are strongly recommended to review the relevant parts of 



II. NON-GAUSSIANITIES IN TWO-FIELD INFLATION 



We are interested in inflation models described by the following action |34j . 135 
S = I d A x 



-g 



±M 2 R - - \e 2b ^g^d,xdvX ~ W{<p, x) 



(1) 



Because of the appearance of b((p), the field x has a non-standard kinetic term. Following 



the notation of slow-roll parameters defined in |2l|, |35 
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p 



W 



~ x 2 V W 



W 



-2b 



Vxx 



e b = 8M 2 b%, 
W 



(2) 



the slow-roll condition can be expressed as C 1, €j, C 1, |^-| 1 with i, j = ip, x- 

As an aside, we mention that model ([1]) is equivalent to the /(x, R) generalized gravity 

(36L 37] when b = —ip/(\/6M p ). But then we find = 4/3, which violates the the slow-roll 

condition. This is a pitfall in treating gen eralized gravity as a two-field model. This pitfall 

can be circumvented by the scheme in [37j. 

Under the slow-roll condition, the background equations of motion are very simple 



3Hip + W tip = 0, 3He 2b X + W t 



0. 



3M 2 H 2 



W. 



Using them one may directly demonstrate 



H 

Jp — e v + e x- 



(3) 



(4) 



Observationally, the most promising probe of primordial non-Gaussianities comes from 
the bispectrum of CMB fluctuations, which is characterized by the non-linear parameter 
/nl(^] , ko, ka). If |/ NL | > 10, it would be detectable by ongoing or planned satellite experi- 
ments I3a, y9l . 
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It has been shown in [24], 34, that the non-linear parameter in two-field inflation 
models can be separated into a momentum dependent term and a momentum independent 
term 



ki,k 2 ,k 3 ) = --/] 



S(ki,k 2 ,k 3 ) 



6 (4) 
r^JNL 



It is also proved in [2 

c(3) 



34J that the first term is always suppressed by the tensor-to-scalar 
ratio, leading to \f^\ <C 1. Hence this term is negligible in observation. For action (pQ), the 
second term 



6 , (4 ) _ N%* N .<p.v. 

g J NL — 



+ 2e 



(K* + e ~ 2b * N2 x*) 



+ e 

2 



(6) 



may be large and deserves a closer look. Here N = f c Hdt is the e-folding number from the 
initial flat hypersurface t = t* to the final comoving hypersurface t = t c . To evaluate we 
will work out the derivatives of N with respect to tp* and x* i n the next section, focusing 
on a class of analytically solvable models. 



III. HUNTING FOR ANALYTICALLY SOLVABLE MODELS 

Making use of equations ([2]), the e-folding number can be cast as 



N 



Hdt 

(W-Q)<p 



dt 



W-Q 



dip — 



Q(<P,X)X 2bM d , 



Qe 



2b 



d\- 



(7) 



Hence Q(<p, x) 1S an arbitrary function of ip and x m principle, because (p/W )(p = e x/W iX = 
— 1/(3H) along any classical trajectory under the slow- roll condition. However, for a given 
W, we have to choose a suitable form of Q so that the integrations defined by Q in f lT9|) 
can be performed. Later on we will fix Q to meet the ansatz ( 12 Up for simplicity. But for the 
moment let us leave it as an arbitrary function of ip and X- It is straightforward to obtain 
the first order partial derivatives 



dN 

dtp* 



dN 

dx* 



W* -Q* W c - Q c dip c 



M%W% 



M*W% dip, 



Q c e 2bc d X c 
W c - Q c dip c 



Q*e 



26, 



dx* 
Q c e 2b 



( W-Q 



M p W *x 



(_Q 



,2b 



dip 

(W-Q\ d X 



dx, 



MIW% d X * 



\M*WJ x d X 
(_Q 



-dip 



dXc 



■x 



dx 

dip,: 



dip 



dip 

\mwj W* 



dx- 



(8) 
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Akin to 24|, |34J, we define an integral of motion C along the trajectory of inflation 



C 



Wf(<p, X ) 
M 2 W tV 



dip 



We 2b f 
M*W~ X 



(9) 



Here the explicit form of f((p, x) is determined by scalar potential W((p, x)- We will give the 
expression of / for some types of potential in this section. If we fix the limits of integration 
to run from t* to t c , then due to the background equations (j3J), 



CI 



Wf 
3M 2 H 



e 2b X 



W„ w, x 



dt = 



(10) 



along classical trajectories under the slow-roll approximation. So the constant C parame- 
terizes the motion off classical trajectories. In order to know dip c /dip*, dx c /d(p*, dip c /dx*, 
dXc/dx* in OS]), we should calculate the first order derivatives of C on the initial flat hyper- 
surface t 



— tlf , 




dC 


w*f* r 
m*w; v + J 


dip* 


dC 


r* ( wf 


dx* 


J \M*W !V 



( Wf 



{ 



M$W >tp 



dx 
dtp* 



dip 



/"( 



( We 2b f \ dcp 
dip, 



M 2 W, X 



-dx, 



■x 



dx 
dx* 



dip 



W*e 2b *f* 
M 2 W* X 



( We 2b f \ dip 



-dx- 



Differentiating with respect to C, it gives 



Wf dip | f f Wf } dx d We 2b f d X f( We 2b f 



M*W tV dC 



\M*W„JdC 



M 2 W :X dC 



/( 



M 2 W, X 



^dx 



(12) 



On large scales, the comoving hypersurface t = t c coincides with the uniform density 
hypersurface. This implies under the slow-roll condition 



W(ip c , X 

whose differentiation with respect to C is 

dip 



const., 



W c 

^ dC 



+ w c 



dXc 
x dC 



0. 



(13) 



(14) 



Combined with (fT2j) on the final comoving surface t = t c , it could give the solution for 
dip c /dC and dx c /dC. This is in general difficult analytically. To overcome the difficulty, we 
introduce an ansatz: 



( Wf 



( We 2b f 
\M 2 W, X 



0. 



(15) 



Although we are free to design the function f(<p,x), the above condition is not always 
satisfiable. We have hunted for analytical models meeting this condition, and found it is 
achievable if W = W(w) with w = U((p) + V(x) or w = U(ip)V(x)- Here W(w), U(ip) and 
V(x) are arbitrary functions of the indicated variables as long as the slow-roll condition is 
satisfied. In this paper, we will pay attention to this situation. But it is never excluded that 
there might be other situations in which dip c /dC and dxc/dC are solvable from (jl2]) and 
( 114"]) . even if ansatz f|T5]) is violated. 
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Ansatz (|T5|) simplifies our discussion significantly. Once it holds, equations (fT2{) and 
lead to 

d<p e 2W C e£e£ 



dC W%ft 



d Xc 2W C e c f c 



while (TTTj) is reduced as 



dc wy< 

dC W*f* 



dp>* M 2 W^ 
dC W*e 2b *f* 
dx~* ~ ~~ MjW* x ' 

As a result, the partial derivatives of iV take the form 

^ + <*-<■ 

fW*-Q*\ fdZ c dQ\ dC lr7C ^ d 2 c 



fl^w* J x \d\* dx*J dtp* 1 ~ d<^*<9x* 



In these equations, we have adopted the notations 



r(w-Q\ d x , ,r(_QfL \ dp. 



dp>* ,<p dC dp* ,x dC dp* ' 

dQ fW c -Q c \ dxcdp c fW*-Q*\ d X * 



dp* V M%W% J xc dC dp* V M*W% J xt dC 

dC r ( W-Q \ (dxY, dC [ c f W-Q \ fx, 
+ dp*J* \m*wJ jXX \dc) V+ dp*]* \M*Wj x d& * 

( Q c e 2b < \ dp c d Xc dC r( Qe 2b \ (dp\\ 

, , , , .. v , 

r«X- 



/ Qe 2b \ dV 



In the above, the expression of Q and its derivatives involve nuisance integrals. To further 
simplify our study, we utilize one more ansatz 



( W-Q \ f_Q 



e 2b 



0. (20) 



In favor of this ansatz, we have Q = and so do its derivatives. 

As was mentioned, ansatz (TT5]) can be satisfied by special forms of potential W(<p, %). 
Now ansatz ( 120]) further constrains the form of W(ip, x) an d °((f). Let us discuss it in 
details case by case. 



A. Case I: W = W(w), w = Ufa) + V( X ) 

For this class of models, according to ffT5"]) . we set 



f = m> (21) 

while condition f[2"0l is met by 

6 = 0, W = Xw a , Q = XVw a - 1 (22) 

or 

*>=-l"U, d ~^ = — l — Q = qe vw W, w . (23) 
2 aw p + qe uw 

Hereafter, as free parameters in our models, A, a, (3, u, p and q are arbitrary real constants. 
The normalization of e 2b is fixed for simplicity. This is always realizable by rescaling the 
field x- 

Taking a = 1, model (I22p recovers the well-studied sum potential 052 ft 

to which 



we will return in subsection IVI Al In subsection IVIDl we will study a specific example of 
non-separable potential that corresponds to a = 2 in (|2"2"|) . 

As will be discussed in subsection IIII C\ there is an equivalence relation between case I 
in this subsection and case II in the next subsection. Models in class I can be transformed 
to those in class II, and vice versa. We will translate model (|23|) to a nicer form (126|) and 
explore it. 

B. Case II: W = W(w), w = U(ip)V( X ) 
For this class of models, we take 

f - ^ (24) 



then condition ( 1151) is satisfied. Condition ( 1201) can be met by 

6 = 0, W = X(\nw) Q , Q = Xilnw^'HnV (25) 

or 

e 2b = U~», ^= ^ v+1 , Q = gw^W, w . (26) 
aw pw + qw v+L 
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We observed that (125} and (j26jl can be obtained from ([2TD, H22) and flU perfectly by 
the following replacement: 

U^kiU, V->kiV, w^lnw. (27) 

In fact, there is a general equivalence relation between case I and case II, on which will be 
elaborated in subsection IIII CI 

Equation ( 126]) dictates W implicitly as a differential equation. To obtain the explicit 
form of W, one should solve the equation. This could be done analytically in some corners 
of the parameter space. For instance, setting v = 0, equation ( )26|) gives 

6 = 0, W = Xw a , Q = (3w a . (28) 

However, if q = 0, it leads to a larger class of model 

W = Xw a } Q = 0, (29) 

leaving b as an arbitrary function of (£. Model ( 12 8 p or ( 12 9 p is separable and can be seen as the 
well-studied product potential j2l|, 34 1. More discussion on models with product potential 
will be given in subsection IVIBI In the case that p = and v ^ 0, we find another model 

e 2b = U a , W = Q = X exp ((3w a ) . (30) 

In subsection IVI D| we will study an example of non-separable potential which corresponds 
to a = 1 in (|30|) . Since v is an arbitrary real constant, equation (126]) can generate many 
other forms of potential W . For example, when p ^ and v = —1, we get a model 

e 2b = U, W = X{w + (3) a , Q = Xf3{w + (3) a -\ (31) 



C. Equivalence between Case I and Case II 

We have classified our models into two categories, corresponding to subsections IIII Al and 
IIII Bl In case I, the potential W(w) is a function of sum w = U(<p) + V{x). In case II, the 
potential W(w) is a function of product w = U(ip)V(x)- After the non-dimensionalization, 
case I can be translated to case II by the transformation 

U^e u , V^e v , w^e w . (32) 

The last relation in (132]) is a corollary of the former ones because UV — > e u+v . On the 
other hand, via transformation (I27p . an arbitrary potential of case I can be transformed to 
that of case II. So the two "cases" are just two different formalisms for studying the same 
models. They are equivalent to each other. We are free to study a model in either formalism 
contingent on the convenience. 

For instance, using the formulae in this section, a model with potential W = Xe~ l3lp x 
and prefactor e 2b = atp 2 can be studied in two different formalisms: 

• Formalism I: W = Xexp(-e w ), b = U/2withw = U+V, U = ln(c^ 2 ), V = ln(/?xV«)- 

• Formalism II: W = Xe~ w , e 2b = U with w = UV, U = ap 2 , V = fix 2 /®- 

But apparently, for this model the calculation will be easier in formalism II. Because the 
dependence of W and b on (p and x is unaltered, the quantization of perturbations is not 
affected by the choice of formalism. For the same reason, the exact dependence of Jnl on 
if and x is the same in both formalisms. 
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IV. MODEL I: W = Xw a , w = U(tp) + V(x), 6 = 



This model is given by f l22|) . which is equivalent to model ( 125]) . Corresponding to this 
model, the number of e-foldings and the integral constant along the inflation trajectory are 



N 
C 



U_ 



M 2 

p 



dip 

IT 



dip + 

dx 
V. v 



v_ 



dx 



(33) 



We have defined the slow- roll parameters in (T5]). In the present case, they are of the form 



a 2 M 2 U% 



a 2 M 2 V 2 
p ,x 



XV, 



XV, 



e b = 0, 



a 2 M 2 U% + V 2 



Xir 



O: 



-M U A M U 





a(a - 


-1) 














o 




a(a - 

Vxx ~ w 2 


-1) 


M 2 V 2 X + 


a r. 

w p ' xx 


Now equations ffTB"]) and 


f|T7|) become 








dtp c 


M 2 U, V V 2 X 




dXc 


M 2 U 2 V v 


dC 


U 2 + V 2 

,<p ^ ,x 


5 

C 


dC ~ ' 


U 2 + V 2 




dC 1 




dC 


1 




dtp, M 2 U 


* ' 

f 


dx* 


M*V*' 



while the function Z defined by (1191) takes the form 

Z 



u 2 v - uv 2 

,<p ,x 



a{U% + V 2 X ) • 

Then we get the partial derivatives of Z c with respect to tp* and x* 

aM p V; x dZ c V2w*A 



aM p U*^ dZ c 



in terms of 



A 



a/2w* dtp* y/2w* dx* 



w c2 e c e c 
^ v x 

w* 2 a 2 e c 



1 + 



a(e% 



r c2 



(34) 



(35) 



(36) 



(37) 



(38) 
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With the above result at hand, it is straightforward to calculate 



N 



,<p* 



w u 



N, 



w v 



aM*V; x : 



N 



aM 2 



x* ~ 

I - ^ ] ,w/ + r+'—A 



n 

2e 



N 



2w* 2 A 



,<p*x* 



aM^V* 



N 



,x*x* 



aM* 



I xx 

2e* 

x 



where for convenience we used notations 

U* + aZ c 



a 
el 



u 



ur 



a A 

(xv + m H A 



V* - aZ c 



to* 



(39) 



(40) 



For these notations, the relation u + v — 1 holds. In the next section, the definitions of u 
and v are different, but the same relation also holds. 

As a result, using formula (jH]) we get the main part of non-linear parameter in this model 



(4) 
K NL 



2 / w v 
- — + — 



+- 



"X 



i _ !^ 

2e* 

f 



au + v 



av + u 



u 



+ - + - a 2 A 



(41) 



The non-linear parameter ff4T|) depends on the exponent a in a complicated manner. 
For the purpose of rough estimation, we assume both u and v are of order unity. This 
assumption is reasonable if U*, V* and w* are of the same order. It is also consistent 
with the relation u + v = 1. Furthermore, motivated by the slow- roll condition and the 
observational constraint on spectral indices, we assume the slow-roll parameters are of order 
0(1O~ 2 ). In saying this we mean all of the slow-roll parameters are of the same order, which 
is a strong but still allowable assumption. After making these assumptions, we can estimate 
the magnitude of (|41l) in three regions according to the value of a. 

Firstly, in the limit a < 1, we have a 2 A ~ ew c2 /w* 2 ~ e. So the third term in curly 
brackets of (f4T|) is of order e -1 u> c2 /u>* 2 , while the other two terms are of order e _1 . Conse- 
quently, we can estimate ~ e/a. It seems that a small value of a could give rise to a 
large non-linear parameter. Specifically, under our assumptions above, if a ~ O(10~ 3 ), then 
the non-linear parameter /nl ~ O(10). However, this limit violates our assumptions. On 
the one hand, we have assumed e v ~ e x ~ r] vx . On the other hand, equations (|34"|) tell us 
e lfi e x /r] 2 x ~ a 2 /(a — l) 2 , which apparently violates our assumption in the limit a <C 1. So 
we cannot use the oversimplified assumptions to estimate the non-linear parameter in this 
limit. 

Secondly, for a > 1, we would have a 2 A ~ aew c2 /w* 2 ~ ae. Then the last term in braces 
of (T4T!) is of order ae~ 1 w c2 /w* 2 . The other terms can be of order ae -1 . After cancelation 
with the prefactor, it leads to the estimation ~ £• That is to say, in this limit, the non- 
linear parameter is independent of a in the leading order and suppressed by the slow-roll 
parameters. 
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The third region is a ~ 0(1). In this region, the non-linear parameter is still suppressed, 

Our conclusion is somewhat unexciting. This model could not generate large non- 
Gaussianities under our simplistic assumptions. However, one should be warned that our 
estimation above relies on two assumptions: u ~ v ~ 0(1) and e ~ f) « 1. Although 
these assumptions are reasonable, they may be avoided in very special circumstances. To 
further look for a large non-Gaussianity with our formula (|4ip . one should give up these 
assumptions and carefully scan the whole parameter space in a consistent way. Generally 
that is an ambitious task if not impossible. But for a specific model of this type, we will 
perform such a scanning in subsection IVI D[ 



V. MODEL II: d\nW/dw = (pw + qw^ 1 )' 1 , w = U(<p)V(x), e 



Jib 



u- 



As we have discussed, model (123]) and model (|26|) are equivalent. Thus it is enough to 
study them in the relatively simpler form, namely in the form fl26|) . For this model, we 
calculated the number of e-foldings and the integral constant along the inflation trajectory 



N 
C 



P U 



m 2 u^ 



qV 



v+1 



M 2 P V X 



U 



M*U, V 



d<p 



V 



M$V X 



d\- 



(42) 



Parallel to section IIVt we also calculated the slow-roll parameters in this model, 
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u v vl 
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Vxx 



2 (p + qw u ) 2 U 2 ' 
2u 2 M 2 U% 



x 2 (p + qw v ) 2 V 2 ' 
M 2 U%V 2 + U U+2 V 2 X 



U 2 



(p + qw v Yw jL 



M 2 [l-p-q(v+l)w v }U. 



M 2 U^ 
(p + qw v )U' 



(p + qw u ) 2 U 2 
M 2 (l - qpw v )U v ' 2 U^V x 

(p + qw u ) 2 w 
M 2 [l - p - q[v + l)w u ]U v V 2 x M 2 U»V XX 



(p + qw u ) 2 V 2 
Subsequently, after obtaining the equations 

dtp c 



+ 



(p + qw u )V 



M 2 UU V V 2 X 



dC 

dXc 
dC 

dC 

d<p* 



U 2 V 2 + U V+2 V 2 

,<p ' ,x 

M 2 P U%VV X 



u%v 2 + U»+ 2 V 2 



u 



M 2 U* ' 

p ,<p 



dC 
dx* 



V* 



M 2 v; x 



(43) 



(44) 
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and 



pU 2 V\ 



u%v 2 



we find by a little computation 



dZ c 



u v+2 v 2 



M p U* u V* x 



dZ c 



V2(p + qw* v )U* d(p* 



V2{p + qw* u )V* dx* 

V2{p + qw* u )A 
M p U* u ' 



Here notation A is different from the one in the previous section, 



A 



U* 2u (p + qw 



cv\2 gC e c 



V X 



[pue 



U c2u (p + qw* u ) 2 e c3 
- c e c - 

~<p x ' ~X'"PV> 



qvw cv e c2 



(4 - 2qvw cu )e c lp e c y + + e> c 



In terms of 



and the relation u + v 



u 



p + Z C U* 



qw 



Z cjj*u 



N 



p + qw* u p + qw* u 

1, once again straightforward calculation gives 

(p + qw* u )U*u (p + qw* u )V*v 



M 2 U* 

p ,f> 



N 



M 2 

p L 



2e* 

*v\2„ 



u 



M 2 U*»V* Y 



A 

pvv H 

e* 



AT 



,V*X* 



2(p + gw^) 2 w*^ 



AT 



x*x* 



^xx 



A 

2e x 



Therefore, the non-linear parameter in this model is 



6,(4) 







2 


















v 2 ' 


+ 




* 




C X L 



+ 



-2 



'XX 



2e* 



u 



2e* 

x 



f + quw u 



+ 



it 



pvv 

V 

el 



A 



Similar to the previous section, we can estimate by assuming u 



(45) 



(46) 



(47) 
(48) 



(49) 



(50) 



0(1) and 



e ~ r\ <C 1. Under these assumptions, the only possibility to generate a large non-linear 
parameter is in the limit v ^> 1. Unfortunately, careful analysis ruled out this possibility. 
Because the assumption e v ~ e& implies pv + qvw v ~ C(l), we find the non-linear parameter 
is not enhanced by v but is suppressed by the slow-roll parameters, ~ e. The same 
suppression applies if v lies in other regions. So we conclude that it is hopeless to generate 
large non-Gaussianities in this model unless one goes beyond the assumptions we made. A 
careful scan of parameter space will be done in subsection IVIEI for a specific model. 
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VI. EXAMPLES 



In sections above, we have generalized the method of 24], |34j and applied it to a larger 
class of models. These models are summarized by equations (|22|) and ff26l) . whose non-linear 
parameters are given by (|4ip and (15771) generally To check our general formulae, we will 
reduce (|41|) and f[5Uj) to previously known limit in subsections IVI Al and IVI Bl The reduced 

HH. In subsections PVTCl IVTDl and IVTEl 



expressions are consistent with the results of [2j 
we will apply our formulae to non-separable examples and scan the full parameter spaces. 

We should stress that all results in this paper are reliable only in the slow-roll region, 
that means at the least e* < 1, < 1, \r]*j \ <C 1 with i, j —if,x- A method free of slow-roll 
condition for some special models has been explored in reference 



A. Additive Potential: W = w, w = U{ip) + V(x), b = 



This potential is obtained from ( 122]) by setting a = 1. The condition b = is necessary 
to guarantee (I2UI) . After taking a = 1, the result in section HVl matches with that in 24 
obviously. 



B. Multiplicative Potential: W = w, w = U{ip)V(x) 

Like equation (129 p . we leave b as an arbitrary function of if, as long as the slow-roll 
parameters (T5]) are small. This is a special limit of section [V] 
Using relations 

-2b 



p=l, q = 0, U u ~- _ 
v = -2 si g n (M si g n f ~jj J Jj-i 



we get the reduced form of non-linear parameter 

6 

5 



r J NL 



+ 



u 2 v 



2e* 




v 4( 






>> X 




\ 2 







(51) 



(52) 



where we have made use of the fact that u + v = 1 as well as the following notations 



M= l_^e 2bc - 2fe *. 



v = JC e 2b c -2b 



(53) 



A 



c c 

e ip € X c 4b c -4b, 
£ c3 



-- sign (6, v ) sign I ) e x 




(54) 



One may compare this formula with [34]]. Note that their definitions of w, v and A are 
slightly different from ours by some factors. Taking these factors into account, the result 
here is in accordance with 34 . 
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C. Non-separable Potential I: W = (oop 2 + fix 2 ) 1 * ', 6 = 



We spend an independent subsection on this model not because of its non-Gaussianity, 
but because it has an elegant relation between the e-folding number and the angle variable 
of fields. For this model, the number of e-foldings from time t during the inflation stage to 
the end of inflation is 

h = (55) 

a(t) v ; AuM 2 AuM 2 y ' 

Note that us can be regarded as sum of squares. Its time derivative gives the Hubble 
parameter ds/dt = —H. So we can follow the standard treatment to parameterize the 
scalars in polar coordinates 

(p = 2M p y/us sin 9, x = 2M p y/us cos 9. (56) 

Rewriting the equations of motion fl3]) in terms of the polar coordinates, we obtain a 
differential relation between s and 9 for the present model, 

. o . dsm 2 9 Rsin 2 9 

sin 2 9 + — — — - = = — (57) 

d\n(vs) Rsin 2 9 + cos 2 9 V ; 



with R— a/ j3. It can be solved out to give 



( s inf?W(«-i) 



N + s e = s = s ^ os ^ 2R/{R _ iy (58) 
At the end of inflation, if the scalars arrive at the bottom of potential, one may simply set 

S e = 0. 

Relation (I5"8"j) is a trivial but useful generalization of Polarski and Starobinsky's relation 



241 l4ll |42| . Recall that Polarski and Starobinsky's relation has been widely used for the 



inflation model with two massive scalar fields, which corresponds to exponent u — 1 in 
the model of this subsection. The simple demonstration above generalized the relation to 
arbitrary u. 

As an application, we evaluate (158]) on the initial flat hypersurface t = t* and then on 
the final comoving hypersurface t = t c , getting the ratio 

Sc _ r^Y^ (^ox r/{r - i] m 

s* vsin^/ \cos9 c J 

which reduces to 

2 / 2 \ R 

r ' / X ' 1 . (60) 



This result can be also achieved from (flU]) directly. 



D. Non-separable Potential II: W = {a<p 2 + f3x 2 ) 2 , 6 = 

Our purpose in this and the next subsections is to examine non-Gaussianities by param- 
eter scanning. Two common assumptions will be used: the e-folding number is fixed to be 
iV = 60 and the inflation is supposed to conclude at the point e° + e£ = 1. 
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Using the latter assumption and the general formulae in section IIVt we find all of the 
relevant quantities can be expressed by e v , e x and R: 



i}tptp - 



(e v + Re x 



XJi 



X 



Vxx 



8M 2 R 2 e x 



(61) 



U 



A 



e * _|_ R e * 
V T Lt X 



+ 



Re* 



(R-l)e c v e c x (e;+Rel 



K + Re c x ) 



(R-l)e^e x (e; + Re 



e^2 



44(4 + Re *x) 2 r 
4(e^ + Re x Y 



s c n + ReZ) 2 ^ 



R 



(62) 



6 f (4) 

g J NL 



/ * 2 1 * 2\ 



-e*V(e* 



.Re! 



+e;e*w(e*u + e*v) + AA{eLu + e*v) 



R 



(63) 



N 



+ xl <pl + xl 



8M 2 

(4 + ^4) 2 



8M p 2 



(64) 



Here we defined R — a/ (3 like the previous subsection. If R — 1, it can be proved that 
-6/4l/5 = (e* + e*)/2 = 1/(2N + 2). Without loss of generality, we will consider the 
parameter region < R < 1. As has been mentioned, from (fTUj) or ( |58|) . one can get relation 
(!60|) . This relation is equivalent to 



R 



efo + Re 



c \ 2(1-R) 
X 



*V T Lt x 



(65) 



If i2 = 1, it gives e^/e; = e x /e* = l/(^ + e*) = iV + 1 and thus r/* x = (2e x + 1)/{2N + 2). 

In the above expressions, there are five parameters: e* , e* e x and i?. The number 
can be reduced by the assumptions we made at the beginning of this section. 3 Firstly, e£ 
and e x can be traded to each other with the relation e£ + ef, = 1. Secondly, since we have 
assumed N = 60, equations (1511) and fl65|) can be used to eliminate two degrees of freedom 
further. Now we see only two parameters are independent, and we choose them to be e c and 
R in the analysis below. The number counting in this way agrees with the fact that (JHJ is a 
first order system under the slow-roll approximation. 



We are very grateful to Christian T. Byrnes for pointing out an error on this issue in an earlier version. 
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As a useful trick, we introduce a dimensionless notation x = xt/xl, then equations (J60 
and (164")) can be reformed as <pl/<Pc = xR anc ^ 



(p 2 x R + x\ x 
SM 2 



N + 



8M 2 



(66) 



Usually the second equation has no analytical expression for the root x, but one may still 
find the root numerically. In the region x > 0, both x and x R increase monotonically from 
zero to infinity, so this equation with respect to x has exactly one positive real root if the 
right hand side is finite. In terms of e£ and R, this equation is of the form 



e c ^x R + R 2 e c x x 



N+ e * + R2e x 



(67) 



Fixing N = 60, the recipe of our numerical simulation is as follows: 



1. Given the values of e c and R in parameter space 0<e°<l,0<i?<l, numerically 



find the root x of equation ( 1671) . where e° = 1 — e c 



2. Compute e* , e*, 77* ^ and according to 



XX 



c \2 



(e c „x R + ReZx) 2 ' 



e x x(e c v + Re x ) 2 
(e c v x R + Re c x x) 2 



(68) 



and equations f lBTj) . 
3. Evaluate —6/^/5 with the formula 
6 



• f (4) 



(4 + ^) 2 



If*" 



[x K + (i2 - l)^ + - (R - l)e 



cl2 



2x R 



3 [x R + {R - l)e 



c l 2 



(i?-l)e 



c|2 



1 - 



(<£ + Re c x ) 2 



R 



x 

-1)6 

x fi + (i2 - 1)4 ifcr - (i2 



+ 



1)4 



-1 2- 



■.(69) 



4. Repeat the above steps to scan the entire parameter space of e c x and R. Due to the 
violation of slow-roll condition, the vicinity of R — should be skipped to avoid 
numerical singularities (see spikes in figure [T]) . 

In a practical simulation, we scan the region < ef, < 1, 0.001 < R < 1 on a uniform 

A 

grid with 101 2 points. Some simulation results are illustrated in figured! When drawing the 
figure, we have imposed the slow-roll condition e* < 0.05, e* < 0.05, r/*^ < 0.05, rf xx < 0.05. 

In the limit R — 1, they are in agreement with the analytical results —6/^/5 = l/(2iV + 2), 
r/* x = (2e x + 1)/(2N + 2). One may also check the results in other limits analytically, such 
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x 1(T 3 



8.7. 
8.6. 
8.5. 




R 



Figure 1: (color online). The non-linear parameter (|69p and slow-roll parameter rf xx (|6ip 
as functions of e c x and R, under the assumptions N = 60 and e° + e° = 1. R is defined 
as R = a/ /3, the ratio of two parameters in the potential of this model. 

as e c — > or e£ — > 1. Theoretically, R — > should correspond to an inflation model driven 
by one field (p. But our method does not apply to that limit, because it would violate the 
slow-roll condition for \. 

From figure [ljwe can see the non-linear parameter is suppressed by slow-roll param- 
eters. Especially, in the neighborhood of R = 0, the spikes of are located at the same 
positions as the spikes of 7/* x . Such a coincidence continues to exist even if one relaxes the 
slow-roll condition. But there is no spike in similar graphs for e* e* and r]* w . Actually, 

these spikes are mainly attributed to the enhancement of /^l anc ^ ^xx ^ 1/-R in the small 
R limit. After the parameter scanning and the numerical simulation, our lesson is that this 
model cannot generate a large non-Gaussianity unless the slow-roll condition breaks down. 
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E. Non-separable Potential III: W = Xe 



-/Vx 2 e 2fc 



This is a special model of ( 126]) with p = 0,z/ = — l,g = — 1. As in the previous subsection, 



we assume N = 60 and e£ + e° 



1. Then from section |V] we get the relations 



2M 2 /3 2 2 
' -<P , 



AM 2 B 



a 

1 hb 

2 V e 



a" 



AT 



In 



C * 

IX 

e c e* 

xv 



%x 



V 



e v e x 



(70) 



c * 
XV 



-4 



*2 c c2 4^20 



2e* 2 

x 



_ f W 



/ * 2 1 * 2\— 2 

(e x M + e^v ) 



e*4M 2 /3 , 

'e y u + e„w + 2e w ) + 2i( £ h - e v) 



a 



For the present model, equation (jTOll gives 



(71) 
(72) 

(73) 

^ - -A (74) 

X X 

If we introduce the notations R = (e?,e*)/(e£,e*), then combining it with equation (j71"j) 



ln(|)=«( X , 



X*) 5 



that is 



In 



c2 * 
X 1 
p*2 
V X 



and the condition e c v + e£ 



1, we can express e* , e* and in terms of e£, e£ and R, 



(R-l)e 



_ 6 t W 

r ■/ NL 



e; = ii%exp 
e* = Re c x exp 
1, /1\1-& C , + i? 2 ^ - l)e c 3 - 2i? 2 (i? - l)Vf 



In , 

5 J " L 2[1 - Re c v + R(R - l)ef} 2 

On the basis of equation (ITUj) . we deduce that ln(l /R)/N should be positive and suppressed 
by slow-roll parameters. In particular, 

* 
* 

^xx 



(75) 



(76) 



* Re % 1 , / 1 
2e; ^lnU 



e x * 



7? 



(77) 
(78) 
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Figure 2: (color online). The non-linear parameter (f76l) and slow- roll parameter 7]* x (|78l) 
as functions of e x and R, under the assumptions N = 60 and e c + e?, = 1. i? is defined 
as i? = (e^ep/(e^ep 5 and it is plotted in logarithmic scale. 

Thus we focus on the region < R < 1. 

As indicated by the above analysis, if we are interested only in the non-linear parameter 
and slow-roll parameters, this model has two free parameters after using our assumptions 
and equations of motion. They will be chosen as e c x and R in our simulation, just like in the 
previous subsection. But we should warn that, compared with the previous subsection, the 
notation R has a distinct meaning in the current subsection. 

The parameter scanning is illustrated by figures |2] and |3j In figure [2J parameter R 
decreases exponentially from 1 to e -60 . In this process, the non-linear parameter grows 
roughly proportional to ln(l/i?) while the slow-roll condition \t]* | <C 1 is violated gradually. 
This phenomenon agrees with equations flTHj) and f!76|) . both of whose amplitude are enhanced 
by the factor ln(l/i?)/iV when R is small. In figure^we find a sharp spike for the non-linear 
parameter in the corner e% — > 0, R — > 1. Figure El is drawn to zoom in this corner, with R 
scaled linearly. As shown by this figure, the spike dwells in a position violating the slow-roll 
condition \rj* J 1. Therefore, the non-linear parameter in this model must be small once 
the slow-roll condition e* C 1, < 1, \r]*A 1 (z, j = (p,x) is imposed. 
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Figure 3: (color online). The non-linear parameter (I76p and slow-roll parameter r\* w 
(|77p as functions of e£ and R near the corner e£ — s- 0, i? — )■ 1, under the assumptions 
N = 60 and e° + e£ = 1. i? is defined as = (e£e*)/(e£,e*), and it is plotted in linear 
scale. In the middle and the lower graphs, the regions with rf w < — 1.5 and rf w < —10.5 
respectively are cut off. 



VII. SUMMARY 



In this paper, we investigated a class of two-field slow-roll inflation models whose non- 
linear parameter is analytically calculable. 

In our convention of notations, we collected some well-known but necessary knowledge in 
section [Til Slightly generalizing the method of 24], |34j , we showed in section IIHI how their 
method could be utilized in a larger class of models satisfying two ansatzes, namely (Tl5*|l 
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and f )20|) . In subsections IIII Al and IIII Bl we proposed models meeting these ansatzes. We 
put our models in the form of W (w ) with w = U ((f) + V(x) in subsection IIII Al and with 
w = U(if)V(x) in subsection IIII Bl At first glance, these are two different classes of models. 
But in fact they are two dual forms of the same class of models, just as proved in subsection 
IIII CI In a succinct form, our models can be summarized by equations (J22J) and (1261) . whose 
non- linear parameters were worked out in sections lIVI and IVl respectively, see equations (jUJ) 
and (1501) . Under simplistic assumptions, we found no large non-Gaussianity in these models. 

As a double check, we reduced the expression (T4T1) for non-linear parameter to the ad- 
ditive potential in subsection IVI A| and fl50|) to multiplicative potential in subsection IVIBI 
The resulting non-linear parameters match with confirming; our calculations. In 

subsection IVI CI for a special class of models, we generalized Polarski and Starobinsky's 
relation f )58|) . For more specific models, we scanned the parameter space to evaluate the 
non- linear parameter, as shown by figures in subsections IVI D l and IVI El In the scanning, we 
assumed the e-folding number = 60 and the inflation terminates at —H/H 2 = 1. For the 
models we studied in subsections IVI Dl and IVI El the non-linear parameter —6/^/5 always 
takes a small positive value under the slow-roll approximation. 
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